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Abstract 
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1 Introduction 


Nonlinearity of field equations ean have pronouneed manifestations in physieal 
predictions from the theory. Examples are ubiquitous and well-known: multi¬ 
ple ground states, static solitons, solitary waves, long lived oscillons, blow-up of 
solutions at a finite time, and many more. Another important effect due to the non¬ 
linearity is that the field can react to the presence of external charges in a rather 
nontrivial way. In consequence, the force between the charges essentially differs 
from naive expectations based on free field models. Very old, yet still interesting 
example of this phenomenon is provided by classical non-Abelian gauge fields 
in the presence of static external charges, see, e.g., [[U, O, (Si. Recently, we 
have studied such this aspect of the nonlinearity using the signum-Gordon model, 
0, 0. The remarkable simplicity of the field equation in this model enables us 
to present exact nonperturbative results, in contradistinction to the case of Yang- 
Mills field. The most striking findings are as follows: there are no Yukawa or 
Coulomb tails of the field in the region far away from charges; the charges are 
totally screened by the field; the force between them exactly vanishes when they 
are separated by a distance that exceeds a certain finite value. In the present paper 
we corroborate these results by solving the two dimensional case. 

Let us briefly remind that the signum-Gordon model involves a single real 
scalar field cp that evolves according to the signum-Gordon equation 

+ g sign 99 = 0 , 

where > 0 is the self-coupling constant. The sign function takes the values 
0, ±1, in particular sign 0 = 0. The Lagrangian has the form 

L = g\^\. 

The field potential U{ip) = 5 f|(^| is V-shaped (| | denotes the modulus). The model 
has been studied in several aspects, such as self-similar solutions (Hi, oscillons 
(or rather breathers) Q, and Q-balls or boson stars in the version with a complex 
scalar field (HI. Let us add that recently the signum-Gordon equation has been 
discussed in a much wider mathematical framework of the theory of partial dif¬ 
ferential equations with compressed solutions (U. In general, the model together 
with its extensions has turned out to be very useful theoretical laboratory in which 
we can study various above mentioned aspects of nonlinear fields. 

In the present paper, similarly as in 0, 0, we regard the signum-Gordon field 
as the mediating field, which generates a force between two point charges coupled 
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to it. In JH, the field generated by the two eharges, as well as the foree, are ex¬ 
actly calculated in one spatial dimension, while in three dimensions approximate 
formulas are given under the assumption that the point charges are close to each 
other. In [Q we study three-body forces in the one dimensional case. Because 
the signum-Gordon equation is nonlinear one, finding its solutions is generally a 
nontrivial task, especially in dimensions larger than one. In the present paper we 
investigate the two dimensional case. Using a formal mathematical connection 
with the planar electrostatics we show that also in this case the point charges can 
be totally screened by a cloud of the scalar field. This screening is easily seen for 
a single charge. The case of two or more charges situated not too far from each 
other is much more difficult. The main problem is the determination of the shape 
of the screening cloud. We have found this shape, as well as an integral formula 
for the scalar field forming the cloud. These results are used in order to derive 
the exact formula for the force between the two identical charges. The force com¬ 
pletely vanishes when the distance between the charges exceeds certain critical 
value (equal to 2i?o(?)> see formula (4) below) that non-analytically depends on 
the strength q of the charges and the self-coupling constant g . 

The plan of our paper is as follows. In Section 2 we establish the presence of 
the total screening in the case of two identical charges in two spatial dimensions. 
We find the exact shape of the screening cloud formed by the scalar field, and we 
obtain the integral formula for the field. Section 3 is devoted to the calculation of 
the force between the charges. Section 4 contains a summary and remarks. 


2 The total screening of charges 


Let us begin from the simple case of a single point-like charge of strength g > 0 
located at the origin in the two dimensional space. The field equation for the static 
field has the form 

A(^ = g sign^j - (1) 

where A is the two dimensional Laplacian. The fundamental solution of the linear 
Poisson equation AG = —q5{x) has the form 


G{x) 


--In , , 

27r In 


\x\ 

^0 


where Iq is a constant. We expect that close to the charge the solution 99 is well 
approximated by G{x), and therefore it is positive. Then the term g sign 99 in (1) 
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has the constant value g, and the exact solution of Eq. (1) has the form 

(p{x) = G{x) + ^ + Co. (2) 

Here cq is another constant. Because it can be included into the constant Iq, we put 
Co = 0. The constant Iq is determined from the requirement that our ^p{x) matches 
the vacuum solution 99 = 0 at a certain radius i?o- The matching conditions have 
the standard form, 

(f = 0, (f' = 0 

on the circle \x\ = Rq. Here ' stands for the derivative with respect to |a;|. The 
final form of solution ( 2 ) reads 




9 1 |T| 0^2 <1 

— In - H —X -, 

27r Ro{q) 4 47r 


( 3 ) 


where 



( 4 ) 


The field (p given by formula (3) has positive values in the circle \x\ < Ro{q), 
hence the assumption that sign^? = +1 is fulfilled. To summarize, the scalar field 
in the presence of the single charge has the form 


Pg{x) 


--2-In 1^1 

27 r Roiq) 


+ fT‘ 


_g_ 

An 


for I x\ < Ro{q), 
for \x\ > Ro{q). 


The field Pq{x) is quite interesting for the following reasons. First, note that its 
dependence on the coupling constant g, as well as on the charge q, is nonanalytic 
at = 0, g = 0. Thus the solution is non perturbative in these parameters. 
Second, the charge q is completely screened by the circular cloud of the field, 
because outside the circle \x\ < Rq the field attains the vacuum value 99 = 0 
exactly. The cloud has the constant charge density g sign^? = g of the opposite 
sign than q. Formula (4) for the radius i?o(Q') is equivalent to the statement that 
the total absolute charge of the circular cloud is equal to q, namely 7 ri?o(g)g = q. 

Now let us turn to less trivial case of two identical point charges of strengths 
q located at the points ±b = (± 6 , 0 ) on the x^ axis in the two dimensional space. 
The distance between charges is given hy d = 2b. Instead of Eq. (1) we now have 

Ap = g sign p — q6{x — b) — q6{x + b). (5) 
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Because of the nonlinear sign {p term, the pertinent solution of this equation is not 
just the sum of appropriately shifted in space solutions (3), unless the distance 
d between the charges exceeds 2Ro{q), in which case the two circular screening 
clouds surrounding the charges do not intersect each other. If the charges are close 
to each other, i.e., d < 2Ro{q), it is not clear whether the total screening is still 
present, and if present, what is the shape of the screening cloud of the field. We 
address these questions in the remaining part of this section. Our findings are 
utilized in the next section, where we calculate the force exerted on the charge 
located atx = b. 

It is clear that when 6 = 0 we again have a single point charge, screened by the 
circular cloud of the field as discussed above, except that now the strength of the 
charge equals 2q instead q. Let us assume for a moment that the total screening 
persists also when 6 > 0, that is that there exists a region S surrounding the two 
charges such that p{x) = 0 for points x lying outside S, or on the boundary (9S, 
and that p{x) > 0 if a; lies inside S. Inside S, the term g sign p = g m Eq. 
(5) provides the constant charge density that screens the two point charges. This 
cloud of charge contributes to p{x) the term 


^doud{x) = 7 ^ / d‘^y In 
2'K Jy, 


The two point charges contribute 


<1 




\x-y\ 
Ro{q) ' 

\x — b\ g j \x + b\ 
Ro{q) 271 Ro{q) ' 


( 6 ) 

( 7 ) 


Note that po{d:) is the solution of (5) when g = 0. The shape of the region S can 
be found from the requirement that we have the total screening, i.e.. 


q^doudix) + po{x) = 0 


( 8 ) 


for all points x outside S. 

As already noticed, Eq. (5) is similar to Poisson equation of ordinary elec¬ 
trostatics. Standard reasoning known from the electrostatics gives the boundary 
conditions for p at (9S: p{x) = 0, dnp{x) = 0 for all x G (9S, where dn denotes 
the derivative in the direction perpendicular to (9S. These conditions are checked 
numerically after we determine the region S. 

The condition (8) is utilized as follows. Eirst, we rewrite the formulas for 
^doudix): ^oix) using the rescaled polar coordinates, introduced as follows: 


X = Ro{q) r 
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sino; 



In particular, r = \x\/RQ{q) and p = |y|/i?o(5')- Thus, 

q /‘2’r rro{a) 2 2p 

‘fcioud{r,0) = — da dpp [21nr + ln(l + ^-^ cos(6'- a))], 

47r^ Jo Jo r r 

Voij', = -- Inr - ^ ln(l + —-cos6') - ^ ln(l + ^ H-cos6'), 

where 

, b 

MqY 

The funetion ro(Q;) gives the radial eoordinate of the boundary of S at the az¬ 
imuthal angle a. The eonsidered set of two point eharges is symmetrie with re- 
speet to refleetions in the both axises, as well as in the origin. We expeet that the 
shape of S refleets these symmetries. Therefore, 


ro(Q;) = ro(27r - a) = ro(7r - a) = ro{n + a). 


Our present task is to determine the funetion ro(Q;). To this end, we first notiee 

that 

ln(l + Y— - - cos(6* — a)) = ln(l — - exp[i(6' — a)]) -|-ln(l — - exp[—*(6^ — a)]). 

lyt ^ ly ly ry 

Expanding the two logarithms on the r.h.s with respeet to p/r, integrating over p, 
and using the symmetry of ro(Q;) we obtain 


^doudij') d) 


47r^ 


J^da rPa) 1„ r - — ^ * 


da rQ^'^{a)e 


ika 


(9) 

This formula is to be eompared with the Fourier series in 9 for (po(r, 9), whieh ean 
be obtained by expanding the logarithms similarly as above, and reads 


= --Inr + 

Itt 


1 + 


k=l 


k 


cos{k9). 


( 10 ) 


Condition (8) is satisfied if 



47r, 


( 11 ) 
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and 


r* 27 r 


da r^+^(a) e*''" = (1 + (-1)"^) {k + 2) n 


( 12 ) 


where k = 1,2,_ Beeause ro(Q;) = ro(7r — a), eonditions (12) are satisfied 

automatieally for odd values of k. For even values of k, we put A; = 2/ in (10) and 
rewrite (11) and (12) as the following set of eonditions for ro(Q;) 


r*27r 


da r, 


2(Z+1) I 
0 


( a ) = 4(/ + l)7rd"'. 


(13) 


^5(«) 


where / = 0,1,2,.... 

Beeause the r.h.s.’s of eonditions (13) depend on we expeet that 
eontains only even powers of d, 

^o('^) ~ /o(q^) + f2{ci) + d'^fi{a) + .... 

The funetions / 2 n(a), n = 0,1, 2,..have Fourier representation of the general 
form 


/2n(a) = y^C2n;fcCOs(2/cQ;), 


(14) 


k=0 


in eomplianee with the symmetries of ro(Q;). The first funetion, /o, gives Tq when 
d^ = 0. Thus, we put = 0 and rg(Q;) = fo{a) in eonditions (13); fo{a) 
has the Fourier form (14). Simple analysis based on the formulas 2 cos(2A;q;) = 
exp(2fA;Q;)+exp(—2f/cQ;) and daexp(ina) = 27r(5„o shows that the eonditions 
for fo{a) are satisfied only if this funetion is eonstant, fo{a) = 2. This result 
is in agreement with the faet that for d = 0 we aetually have the single point 
eharge of the strength 2q, whieh has the eireular sereening eloud of the radius 
i?o(2g) = 2i?g(g), i.e., = 2 for our resealed radial eoordinate. 

In order to determine the funetion f 2 (a), we differentiate eonditions (13) with 
respeet to d^, and we put d^ = 0, fo = 2. This gives the eonditions 


r*27r 


/ dtt e^*'"/2(a) = 27r dn, 

do 

where I = 0,1, 2,.... It follows that /2 = 2 cos(2q;). Taking higher derivatives 
of (13) with respeet to d? and performing similar ealeulations as above, we have 
found that /4 = 0, /e = 0, /s = 0. With sueh partial results, we have made the 
edueated guess that / 2 „ = 0 for all n>2, i.e., that 


rg(Q;) = 2 [1 + d^ cos(2q;) . 


(15) 
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It turns out that indeed, rl{a) given by this formula satisfies all conditions (13). 
The pertinent integration on the l.h.s. of (13) is elementary. 

The form of formula (15) implies that it holds only if d < 1, because Tq > 0. 
In the case d = 1, i.e., |6| = Ro{q), this formula gives two circles which touch 
each other at the origin, and have the point charges at their centers. For d > 1 
each charge has its own circular screening cloud, separated from the other. 

It remains to check whether all values of the total field (^(x) = ^Pdoudix) + 
ifoix) are strictly positive inside the region S, as it has been assumed. We can do 
this only numerically. It is a rather straightforward computation since the integral 
(6) giving (f cloud has the already known compact domain S. Moreover, the two 
dimensional integral can be reduced to one-dimensional one over the boundary 
of S, see formula (19) below. The numerical results corroborate our assumption. 
They also show that cp and are continuous on OS. 

Thus, we conclude that the pertinent solution of Eq. (5), denoted below as 
p>qq{x), has the following form 


P>qq{.x) 


Vdoudix) + (po{x) for f e S, 

0 for X ^ S, 


(16) 


where the compact region S has the boundary with the rescaled radial coordinate 
ro(a!) given by formula (15). Example plots of levels of (Pqq{x) are presented in 
Figs. 1 and 2. 


3 The force exerted on the charge located aix = b 


The method of calculating the force exerted on an external charge coupled to a 
field is discussed in detail in dUl. Adapting it to the case at hand, we extract from 
the field (pqq the logarithmically divergent proper field of the charge located at 

X = b, 

Q , \x — ^ 


Here 


m(x) 


g , lx + 61 

— In- 

2vr Ro{q) 


3 “ ^doud{x') 
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Figure 1: The contour plot of (pqq ford = 1/2. The outermost contour corresponds to 
if = 0 and it is given by ro(Q!), formula (15). This contour encircles the compact domain 
S. The horizontal and vertical axises correspond to /Ro{q), x‘^/Ro{q), respectively. 



Figure 2: The contour plot of ipqq for d = 3/4. The meaning of the lines is the same as 
in Fig. 1. The picture shows the deformation of outer contours in vicinity of the vertical 
line = 0, that ultimately leads to the breakup of the cloud into two non overlapping 
circular clouds when d= 1. 
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is smooth at the point x = b. Next, we ealeulate the time rate of the transfer of 
momentum from the field to the eharge loeated atx = b. It is equal to the foree F 
exerted on that eharge, see dH for details and examples. It turns out that 

F=qVu\^^^, (17) 

where Vm = {du/dx^, du/dx^). We see that the proper field of the eharge does 
not eontribute to the foree exerted on that eharge, as expeeted. The symmetry 
of the set of eharges implies that = 0. The formula for the non vanishing 
eomponent F^ ean be rewritten in the polar eoordinates as 


du 

q du 

^ 01-1 
o\x\ 

|lr=b ^o(g) dr 


(18) 


Computation of the foree from formula (18) is signifieantly simplified when 
we rewrite q^doud inside the region E as a line integral over the boundary OS. First, 
we write i-pcioud in the form 


^cloudi,^') 


d‘^y {dyiW'^ 


dy.W^] 


where 


W\x,y) = ^^{x^ -y^){hi 


k^n^ {S-yf 


Rl{q) 


-1), 


dy.W^-dy.W^ = 2\n ^ 

Rl{q) 

and is the antisymmetrie symbol, i,k = 1,2. Next we apply Stokes theorem. 
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Sir 


9 ^cioud{x) = — (p da d^yoia) W{x,yo{a)), 


dT. 


(19) 


The eurve yo{a), defined by 


yo{a) = Ro{q) ro{a) 


cos a 
sin a 


represents OS. Using the eontour representation (19) we obtain 


cloud 


Or 


0=0 

r=d 


gRliq) 

Att 


r-27r 


l^2^T 


da hi{a) h 2 {a) + 


da / 13 (a) / 14 (a) 


, ( 20 ) 
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where 


hi{a) = 


d — ro(a) cos a 


— 2d ro(a) cos a + rl{a) ’ 


/i 2 («) = rl{a) — d ro(a) cos a + 


2d^sma sin(2a) 
ro(a) 


h^i^a) = \n[d? — 2d ro(a) cos a + rl{a)] — 1, 
2 sinasin(2a) 


hdia) = 


ro{a) 


ro{a) cos a, 


ro(a) is given by formula (15), 0 < (i < 1. 

In spite of their appearanee, the integrals in (20) are elementary. First, no- 
tiee that h 4 {a) = —da{ro{a) sin a). Therefore, in the seeond integral in (20) we 
may use integration by parts in order to eliminate the logarithm funetion. See¬ 
ond, we write the sin and cos functions in terms of exp(±za) and observe that 
da exp(ika) = 0 for any integer k except fc = 0. This eliminates many 
terms, in particular all terms proportional to \Jl + cos(2a). This last function 
is thought of as the Taylor series with respect to cos(2a), obviously it contains 
only even powers of exp(±zQ!). Finally, we use the integral 


f-2-w 


da 


Iq 4 + # + 4^2 cos(2a) 


27r 


and formula (4) for Ro{q). We obtain 


cloud 


dr 


= d{2 — d^ 

dvr 


( 21 ) 


The non vanishing component of the force exerted on the charge located at the 
point = 6, = 0 is computed from formula (18), 


= 


Airb 


1 - 






Airb 


1 - 


7rgb^ 


( 22 ) 


Here we have returned to the original coordinate b = RQ{q)d. The force is at¬ 
tractive one, as expected from the scalar field. Formula (22) holds for b < Ro{q), 
outside this range of b the force vanishes 0. 

' Rq (q) vanishes in the limit q ^ 0. Therefore, the force also vanishes in this limit, in spite of 
the fact that the r.h.s. of formula (22) gives (A. The point is that in that limit formula (22) 

is not valid for any b — the force vanishes for all b > Ro{q). 
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The factor —q^/Anh in formula (22) represents the standard two dimensional 
Coulomb force characteristic for the free field {g = 0). This force dominates at 
short distances, I -C Iq, also in the case of self-interacting field. Here Z = 26 is 
the distance between the charges, and Zq = 2 i?o(?) is the critical distance between 
them. At distances Z ^ Zq the force becomes weaker, and it exactly vanishes for 
Z > Zq. When the distance Z reaches Zq, the screening cloud splits into two non¬ 
overlapping circular clouds that screen each charge separately. 

Formula (22) for the force can easily be generalized to the case one particle is 
located at the point x and the other one at y. It is sufficient to substitute b = 1/2, 
where Z = \x — y\, and to include the unit vector n = {x — y)/\x — y\ directed from 
y to X. Thus, the force exerted on the particle located at x is given by formula 



(23) 


^1 < Zo, otherwise F{x) = 0. This force possesses the potential Uqq 
(f) = -dsUqq, namely 


if Z = IT — ^1 ■ 
such that F{x) 



if Z < Zo, otherwise Uqq = 0. 

4 Summary and remarks 

1. As the most valuable result of our work we regard the exact formula (23) for 
the force exerted by one point charge on the other. This force vanishes quadrati- 
cally when the distance Z between the charges approaches (from below) the critical 
value Zo = 2^q/'xg. In the one dimensional case investigated in Q, the force van¬ 
ishes linearly, F = {1 — a/a^)/2, where a is the distance between the charges, 

and the critical distance a* is given by the formula a* = q/g (we keep the same 
form of the field equation in all dimensions). The situation in three dimensional 
case remains to be investigated because in [|4|| we have been able to compute the 
force only approximately, under the assumption that the charges are close to each 
other. 

The field g) is the sum of g)cioud, given by the exact integral formula (19), and 
of the two logarithmic Coulomb terms. In the one dimensional case the pertinent 
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field is known analytically BUl, while in three dimensions we only have an approx¬ 
imate formula for it. Note that formula (19) provides quite convenient starting 
point for numerical computations of the field. 

2. The total screening of the charges coupled to the signum-Gordon scalar 
field may resemble the phenomenon of total screening of external color charges 
interacting with a classical Yang-Mills field [|T]|, [13. One should however note 
that there are several differences: in the Yang-Mills case the total screening is 
proven for external charges that are spatially extended; the screening field is time- 
dependent; the complete screening appears only in a certain limit. Moreover, the 
analytic form of the screening field is not known. It is not clear to us whether the 
two cases of the total screening are interrelated in some way. 

3. In paper [[TOl it is shown that a classical scalar field in the presence of point¬ 
like external charges can be utilized in order to unravel certain essential features 
of the corresponding quantum theory, like asymptotic freedom or triviality, de¬ 
pending on the sign of coupling constant. The main role is played by a classical 
perturbative solution that has the form of a formal series in powers of the cou¬ 
pling constant. The field studied there is the real, massless scalar field with the 
self-interaction of the form Similar investigations in the cases of a massive 
scalar field and Yang-Mills field are presented in ifTTl . [fT3 . respectively. We think 
it would be interesting to employ such a method to the signum-Gordon model, in 
particular because very little is known about the properties of the quantum version 
of this model. 
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